Josephson current through a long quantum wire 
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The dc Josephson current through a long SNS junction receives contributions from both Andreev 
bound states localized in the normal region as well as from scattering states incoming from the 
superconducting leads. We show that in the limit of a long junction, this current, at low temper- 
atures, can be expressed entirely in terms of properties of the Andreev bound states at the Fermi 
energy: the normal and Andreev reflection amplitudes. This has important implications for treating 
interactions in such systems. 



As was shown by Josephson [l| a current can pass be- 
tween two superconductors separated by a normal mate- 
rial, even with zero potential difference. At temperature 
T, this Josephson current is determined by the difference 
of the phase of the order parameter in the two supercon- 
ductors, X- 



ax 



(1) 



where F is the free energy [2]. Using the Bardcen- 
Cooper-SchricfFcr (BCS) approximation in the supercon- 
ducting leads and ignoring interactions in the normal re- 
gion, I(x) can be expressed as a sum over single quasi- 
particle energy levels, E n : 



I[ X ;T] = 2eJ2f(E r , 



dE n 
dx 



(2) 



where f(E) = \ j\e E l T + 1] is the Fermi function and we 
measure energies from the chemical potential. In general, 
these states are of at least two distinct forms. There are 
Andreev bound states (ABS) [3J, with energies \E\ < A, 
where A is the superconducting gap, which are local- 
ized in the normal region and whose wavefunctions de- 
cay exponentially into the superconducting leads. There 
are also scattering states (SS), with energies \E\ > A 
corresponding to waves coming in from infinity in the 
superconducting leads and being reflected and transmit- 
ted. If the bottom of the band in the normal material is 
lower than the bottom of the band in the superconducting 
leads, there are, in addition, normal bound states, local- 
ized in the normal region which also decay exponentially 
in the leads. In general, all types of states contribute to 
the Josephson current. 

The limit of a long narrow normal region was consid- 
ered in [ij , using a nearest neighbor tight-binding model 
and initially ignoring interactions. In particular, it was 
assumed that the length of the normal region, I, was 
much greater than the coherence length, or equivalently 
than the finite size gap, ttvf/£ *C A where Vf is the 
Fermi velocity in the normal region. Furthermore, only 
T = was considered. In this limit it appears natural to 



integrate out the gapped superconductors and derive an 
effective Hamiltonian for the normal region, with local 
pairing interactions induced by the proximity effect at 
its boundaries. Such as effective Hamiltonian was used 
to derive the Josephson current. In this approach, only 
ABS's are considered. Due to a remarkable cancellation 
between pairs of ABS's it was found that the current, 
to order l/£, could be expressed in terms of scattering 
amplitudes at the Fermi energy only. 

This approach was called into question by the results of 
There it was verified that the ABS's gave the entire 
Josephson current for long junctions in the unrealistic 
limit A > 2J where 4J is the bandwidth in the nor- 
mal region. However, numerical results for intermediate 
length junctions seemed to suggest a significant contri- 
bution from SS's for A/(2J) < 1. 

In this letter we study general models of long non- 
interacting SNS junctions without integrating out the 
superconducting leads. We prove that, for vpjt and 
T <C A, the Josephson current can indeed be expressed 
in terms of data at the Fermi level only. We emphasize 
that states far from the Fermi energy make large contri- 
butions to the current; it is just that these nearly cancel 
for large £, leading to our result. 

This finding is important because integrating out the 
superconductors provides a powerful method for includ- 
ing interaction effects in the normal region, based on 
boundary conformal field theory techniques [4j. (See also 
@.) While [1, IU only considered the dc Josephson cur- 
rent, the techniques introduced there can be extended 
to the ac case by allowing for the phase of the boundary 
pairing interactions to evolve linearly in time, x = eVi, 
where V is the voltage difference. A possible experimen- 
tal realization of such a long SNS junction might be pro- 
vided by a carbon nanotube between bulk superconduc- 
tors. Using vp » 8.1 x 10 5 m/s, irvp/l .5 meV for £= 
3 microns. Thus, obtaining sufficiently long clean nan- 
otubes coupled to sufficiently high T c superconductors to 
satisfy itvf/£ *C A may be near the limits of current 
nanotechnology. 

Our proof is based on expressing the total current, 
summed over all types of states, as a contour integral in 
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the complex energy plane, involving the S-matrix. In par- 
ticular, we derive this contour integral for the T = case 
and observe, using analytic properties of the S-matrix, 
that the contour can be deformed into an integral along 
the entire imaginary S-axis. We then show that, in the 
large i limit, Vpjl, <C A, the integral is determined by 
the S-matrix at the Fermi energy, and apply our gen- 
eral result to the "Blonder- Tinkham-Klapwijk (BTK) 
model" , obtaining an explicit formula for the current for 
vf/£ "C A. We extend our contour methods to finite 
T, by expressing the resulting current in terms of a sum 
along the imaginary energy axis at the Matsubara fre- 
quencies, E = iu> n = i2ir(n + 1/2)T. We find that the 
current vanishes exponentially when T 3> vf/L Finally, 
we mention that our approach may readily be extended 
to tight-binding models, such as the one discussed in 
whose results for the current we recover when T — and 
v F /£ < A. 

We consider a general SNS model in the non- 
interacting, BCS approximation with a gap function 
A(ie) of magnitude A at |x| — » oo and in the central 
region. We also include a normal potential, V(x) which 
vanishes at \x\ —¥ oo. A 4x4 transmission matrix, M may 
be defined which relates the asymptotic wave-function in 
the S regions at x — > ±oo, A + = MA" with Bogoliubov- 
DeGennes wave-function obeying (for (x — > ±oo), respec- 
tively) : 



u(x) 
v{x) 



cos(*/2) 
-e***/ 2 sin(tf/2) 



|± e i/3 p x 



cos(*/2) 



sin(*/2) 



[Afe-^ hX 



+ Afe lhx ] , 



A ± e -i/3„x] 



(3) 



where the particle and hole momenta, for energy E, are 
Pp/h = ^ m s{^ ± {E 2 — A 2 )2} ; with ms the electron 
effective mass in the S regions and /j, the chemical po- 
tential, and ^ = — arcsin(A/_E). (For simplicity, we 
take an energy-independent gap, A, but our results can 
be extended to more realistic models.) Eq. ^ applies 
also to the ABS regime, in which E 2 — A 2 < 0. In this 
case, the phases of the arguments of the complex square 
root functions are always chosen so that Im(/3 P ) > and 
Im{(3h) < 0). The S-matrix, which expresses outgoing 
waves (Af , A$ , A% , Aj) in terms of incoming waves can 
be expressed in terms of M. In particular: 

detLSl = M M M 3,3 - M h3 M 3A F{E ]X ) 

M 2 . 2 M AA -M 2A M^ 2 G(E;x) ' 

In Eq. ((4]), F(E;x) and G{E;x) must be regarded as 
functions of E in the complex E'-plane. We chose them 
to obey several convenient properties, which are crucial 
for our derivation: 

i) They are always finite for finite E. This can be easily 
achieved by shifting poles of Q into zeroes of T and vice 
versa; 

ii) They have no common zeroes. [Possible common 



zeroes (e.g. Eq), could always be cancelled by a re- 
definition: F{E- X ) HE;x)/(E - Eo), g(E; X ) -> 
g(E; X )/(E - E ), without changing Eq. ©]. 
Hi) J-(E;x) = G*(E;x)- Here this equation refers to 
complex conjugating the function without complex conju- 
gating its argument, E. This condition is consistent with 
the requirement that |det[5]| = 1 for scattering states. 

iv) Q(E; x) can be denned to have branch cuts along the 
real £7-axis, corresponding to the nonzero density of scat- 
tering states in the leads. This is due to the fact that 
Q (E; x) depends on E via f3 p and fih and that they be- 
come double- valued functions of E, for \E\ > A. 

v) d x \a.Q{E-,x) vanishes rapidily at \E\ — ) oo along any 
ray not parallel to the real axis. This condition is crucial 
to allow for pertinently deforming the integration path 
in the energy plane, when computing /^[x]- 

vi) Q(E; x) is rea l m the bound state region: the real axis 
with -A < E < A. 

These conditions appear to determine F{E; x) and 
Q (E; x) uniquely except for an overall multiplicative con- 
stant factor. Moreover, they imply that zeroes of G{E; x) 
correspond to poles of det[5]. These conditions imply 
that there are no poles of dct[5] off the real axis. We 
are actually dealing with a 2-sheeted Riemann surface, 
due to the branch cuts. We may regard J-(E; x) as be- 
ing G(E; x) on the second sheet of the Riemann surface. 
With the definition we gave of j3 p (E), fih{E), the zeroes 
of G(E;x), corresponding to poles of detfS 1 ], occur ei- 
ther on the real axis, or else off-axis on the second sheet 
of the Riemann surface. This property of det[5] follows 
from general principles. Since the S'-matrix can be de- 
rived from the retarded Green's function it should have 
no singularities in the upper half plane. Since 



Q{E-x)=G*{E;x) 



(5) 



in the BS region, we can use the Schwartz reflection prin- 
ciple to define its unique analytic continuation to the en- 
tire first sheet of the Riemann surface, where it obeys Eq. 
(0. Thus if Q(E;x) had a zero in the lower half-plane, 
at -Eo; it would have to have a twin at energy Eq in the 
upper half-plane. This would violate this basic property 
of S telling us that no such zeroes exist. 

The ABS's correspond to poles of the ^-matrix and 
therefore to the zeroes of Q(E). This allows us to write 
the contribution to the ground state energy from ABS's 
as: 



hj 



(0) 
ABS 



dElnG(£;x) 



(6) 



where the contour Tabs i n the complex energy plane 
surrounds the negative e nergy AB S's. In the SS energy 
regime, for A < \E\ < \J + A 2 (SSI energy regime), 
the S'-matrix is unitary and J- = G* . det S(E) can be 
written in terms of the particle and hole phase shifts for 
the 2 scattering solutions at each energy, (Tp,a^, respec- 
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tively: 



a) 



Im(E) 



detS=n«P[i(o?-^)] 

a=l 



(7) 



On the other hand, for \E\ > ^fj, 2 + A 2 (SS2 energy 
regime), only particles propagate, but we find Eq. ([7]) 
remains true with cr^ set to zero. In this way, we obtain 
the total contribution to the ground state energy from 
scattering states in the form: 



E { Q } = e° QQ — I dEln 



1 



J ss 



z ss 



2ni 



G*(E;x) 



(8) 



where e® ss is independent of \- Using Eq. ([5]), the second 
term in Eq. ((SJ) can be also be written as a contour inte- 
gral in the complex energy plane, like Eq. ((5]), with the 
contour now running on both sides of the branch cut in Q 
along the real E-axis from — oo to — A. These two terms 
can be combined, allowing us to write a simple unified 
formula for the zero-temperature Josephson current: 



/ (0) M = - 



2e 
27ri 



dEd x {\nG(E; X )} 



(9) 



where the contour T runs infinitcsimally above and below 
the negative E axis. Due to the properties of Q(E;x) 
discussed above, T can be deformed into a single line 
running along the imaginary .E-axis from — oo to oo. [See 
Fig. ©•] 

We now assume that our system is made of two su- 
perconductors at phase difference x separated by a long 
central normal region C of length I as sketched in Fig. 
([2]). We assume that A(x) makes an abrupt transition 
from Ae ±4X / 2 to in the central region and V{x) makes 
an abrupt transition from in the leads to Vc in the cen- 
tral region. Here "abrupt" means rapid on the scale of I. 
We assume Vc < M so that the central region is metal- 
lic. The mass in C is written as m. The wave-functions 
in the central region, far from the interfaces, may be 
written u(x) = Ci exp(ia p a;) + C2 exp(— ia p x), v(x) = 
C3 exp(— iahx)+C4 exp(ia/ 1 a;) with a p /h = {2m(ii— Vc± 
E)}^ . We may define transmission matrices for the left 
and right interfaces, L and R by C = LA~ , A + = RC in 
terms of which M = R-M c -L, where M c is the transmis- 
sion matrix of C, given by the diagonal matrix with en- 
tries [exp(za p £), exp(— ia v t), exp(— iotht), exp(za? l £)]. For 
large £, the integral along the imaginary i?-axis is dom- 
inated by energies of 0(vf/£) and the leading term in 
l/£ is determined by the L and R transmission matri- 
ces at E = 0. These in turn can be expressed in terms 
of the normal and Andreev reflection amplitudes at zero 
energy for particles and holes at the left-hand (L) and 



at the right-hand interface (R), N' 
result, one obtains 



p/h 

L/R 



and A 



p/h 

L/R- 



AevF 



As a 



(10) 
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FIG. 1. Sketch of the deformation of the integration path F 
used in Eq.© to get to Eq.JlO]): 

a) The integration path T sketched as F = Tabs U Tss , with 
Tabs running around poles corresponding to ABS's and Tss 
going around the energy interval corresponding to negative- 
energy SS's; 

b) Adding the arcs E + , E~ , whose contribution to the integral 
is zero as their radius is sent to 00 (see text for the discussion), 
allows for trading the integral over T for an integral over the 
imaginary axis. 
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FIG. 2. Sketch V(x),A{x) in the SNS system. The L- and 
R-interfaces are assumed to be "sharp", that is, V(x),A(x) 
and the single-particle mass are assumed to vary over typical 
length scales <C I- In this case scattering is localized at the 
interfaces and is fully encoded within the particle- and hole- 
normal- and Andreev- scattering amplitudes at the interfaces, 



L/W L/R' 



with d(x) = &iccos{Re[N p R N p L e 2iaFt + A p R A h L e lx ]}, and 
ap being the Fermi momentum. This leads to Ishi's 
sawtooth current [9] in the limit of perfect Andreev scat- 
tering at zero energy (N P , R — 0, \A P L ^ l R \ — 1) and a 
complete suppression of the Josephson current when the 
Andreev reflection amplitude vanishes at either interface. 
As an example, we explicitly compute 1^ [x] in a 
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model system whose S-N interfaces may be thought of 
as a generalization of the one studied in [Io| . that is, we 
assume that m is uniform, A(x) abruptly changes at the 
S-N interfaces and 

V{x) = Vo[5(x) + 5{x - £)} + V c 6(x)6(£ - x) (11) 

where 8 and 9 are the Dirac delta-function and Heavyside 
step function respectively 



11]. Eq. 



i?(x) = arccos^Re 
(/3f — otp — iZ)(/3p + aj? + iZ 



-2cos(x) 



TOj) now gives: 



F + \ 
2 



-iZ 



4 + \fo - iZ\ 



„2ia F i 



(12) 

where qj? = ^/2ra{\i — Vc), Pf — {2m(/i + iA)}i and 
Z = 2mVo. Note that two conditions must be satisfied for 
perfect Andreev reflection and hence a sawtooth current: 

Z = lmP F = ^mtvV + A 2 " - n), V c = -Z 2 /(2m). 
This is actually consistent with the results obtained in 
Ref.[J] within a tight-binding approach to the calcula- 
tion of the dc Josephson current across an SNS sys- 
tem: getting perfect Andreev reflection requires, in gen- 
eral, fine-tuning two system parameters. On the other 
hand, allowing the central region to host an effectively 
attractive interaction between electrons should allow for 
the system to dynamically self-tune to theperfect An- 
dreev reflection point, as £ becomes large |J|. Our gen- 
eral result, Eq. flTUl) . also holds for tight-binding mod- 
els where the spectrum may also include normal bound 
states. Indeed, for the nearest neighbor tight-binding 
model at half- filling studied in jj], one obtains i?(x) = 

r 4A|co S ( x ) + (l-A 2 B ) 2 " 

arccos ^ (i+A%)* 

coupling Ab defined as in [4jj. 

Eg. ((TO)) gives the main formula for the dc Joseph- 
son current in Eq.(3.13) of derived using a low energy 
effective boundary mode, taken in the symmetric case 
(A L = A R ). 

It is easy to generalize our contour result to compute 
the dc Josephson current at finite temperature T, T\. 
Now the contour T gets deformed into a sum of circles 
around the poles of the Fermi function at uj n — 2irT(n + 
1/2), yielding 

i[x;T] 



with the effective boundary 
Using this formula for 



2eT £ 

= — OC 

d x Re[A R A h L e^\ 



71— — OG 



cosh ( 2fiJai 



Re[N R Nle 2ia ^ e + A R A h L e^} 



(13) 



Of course, this gives our T — result at T <C vf/£ 
where we may approximate the sum by an integral. At 
T 3> vp ft we may approximate the sum by the two terms 
with uj n — ±7rT: 



I[ X ;T] » %eTe~^ Tll l VF d x Re[A p R A h L e lx }+0 



-6-rrTl/vF 



This becomes exponentially small when T 3> vf/£- 

By using analytic properties of the scattering matrix 
for an SNS system, we have expressed the Josephson cur- 
rent, which has contributions from both bound and scat- 
tering states, as a single contour integral in the complex 
energy plane. In the limit of a long central region, we 
then proved that the current can be expressed in terms 
of properties of the Andreev bound states at the Fermi 
energy only: namely the normal and Andreev scatter- 
ing amplitudes. The result holds at finite temperature 
provided that T,vf/( <C A. This result shows that the 
Josephson current is a universal quantity, in this long 
length low temperature limit, and justifies the low en- 
ergy Hamiltonian approach used in which was cru- 
cial for treating Luttinger liquid interaction effects. It 
also paves the way towards an extension to the non- 
equilibrium (AC) Josephson effects. 

Explicit derivation of Eq.(9) and of the 
formula for $(%) 

The starting point to derive Eq.(9) and, specifically, 
the explicit formula for $(x) in terms of the normal- and 
of the Andreev-reflection amplitudes at the Fermi level, 
is writing detfS 1 ] (that is, J-(E; x) and Q(E;x)) m terms 
of the R- and L-matrix elements. By definition of T and 
Q and of the R- and L-transmission matrices, one finds 
that the following explicit formulas hold: 



T (E; X ) = F ,o(x; E) + F 1A (E)S a "- a ^ 

+ F_i - 1 ( L E)e- i[a »- ah]e + F 1 - 1 (E)e ila * +ah]e 

+ F^ l {E)e- i[ap+ah]e 

g (E; X ) = Go,o(x; E) + G 1 , 1 (E)e i ^-^ e 

+ G-i^i{E)e- i[a '~ ahV + G 1 ,- 1 (E)e i[a » +ah]e 

+ G- 1A (E)e- l[a ? +ah]e , (15) 

with 



FoAx) 


= —(Ri,iRa,2 


— i?i ) 2-R3,l)(-£/l 1 3-^2,1 


— L\. 1-^2,3) 




— (Ri, 3R3A — 


^1,4-^3,3) (-^3,3-^4,1 — 


^3,1-^4,3) 


Fi,i 


= (Ri,iRs,3 — 


-Rl,3-R3,l)(-C'l,l-^3,3 — 


^1,3^3,1) 


F-h-i 


= (Ri, 2R34 — 


^1,4-^3,2) (-^2,1-^4,3 — 


^4,1-^2,3) 


Fi,-i 


= (Ri,iRs,4 — 


-R3,l-Rl,4)(-C'l,l-^4,3 — 


^1,3^4,1) 


F-i,i 


= (Ri, 2^3,3 — 


-R3,2-Rl : 3)(^2,li3,3 — 


^2,3^3,1016) 



(14) 



and 



GqAx) = — (-^2, 1-^4,2 — ^2,2-^4,1) (-^1,4-^2,2 — ii, 2-^2,4) 

— {R2, 3-^4,4 — ^2,4^4,3X^3, 4^4, 2 — ^3,2^4,4) 
Gi.i — (R2, 1-R4.3 — ^2,3-^4,1) (-^1,2-^3,4 — ^3,2-^1,4) 
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G-i-i = (i?2,2-R4,4 — -^2, 4^*4,2) {L2, 2L4, 4 — L2, 4-^4,2) 
= (-R2, 1-^4,4 — -^4,1-^2,4) (-^1,2-^4,4 — il,4-^4,2) 
= (i?2, 2-^4,3 — -^4,2-^2,3) (-^2, 2-^3,4 — ^3, 2-^2,43- 7) 

and the explicit dependence upon x set according to the 
definition of the i?- and L-matrices. Using general prop- 
erties of the transmission matrices, arising from the con- 
tinuity equation for probability current, it is not difficult 
to show that Eqs. (|16ll7j) are consistent with the identity 
F{E;x) = G*(E;x), that G{E; X ) is real for real E and 
—A < E < A, and that the branch cuts of Q(E;x) ne 
on the real axis, from E — > — oo to i? = - A and from 
E = A to E — > 00. From Eqs. (fl~5|) . one then sees that 



d x \nG(E; X ) 



d x G , (E;x) 

Q{E;x) 



(18) 



Because, as E goes to infinity along a radius that does 
not overlap with the real axis, either the imaginary part 
of a p , or the imaginary part of ah, goes to — oo, from 
Eq.(fl~8"|) we find that d x In Q(E; x) exponentially vanishes 
as \E\ — > oo off the real axis. Using Eq.(8) for the current 
and taking into account that, for large t, contributions 
to the integral with \w\ > V are strongly suppressed, we 
may approximate a p and ah as 



a p « aF + iw 



ah ~ ap — ilj 



2(/x - V C ) 



2(/x - V C ) 



(19) 



with ap — y/2m(n — Vcj. Using Eqs. ([T9l) . the zero- 
temperature Josephson current to leading order in l~ x 
is then given by 



+ G_ 1 ^ 1 e- 2M ^ + G ,o(x)}" 1 



2ia F £ 



(20) 



with the coefficients G ai t being defined as the coefficients 
G at b evaluated at lo = 0, that is, setting a p = ah = aF, 
ftp = (3^ — {2ms[n + iA]}5. Computing the integral in 
Eq.(f2"U|). one eventually finds out 



In order to prove that Eq. (|2"Tj) yields Eq.(9) of the main 
text, we have to rewrite Eq.fl22p in terms of the normal- 
and Andreev-scattering amplitudes at the Fermi level. 
To do so, we relate the scattering amplitudes at both 
interfaces to the R- and L-matrix elements. This may be 
readily done starting from the definition of the normal- 
and Andreev-scattering amplitudes. The result is 



A p R {E) 
n r(E) 
A h R (E) 



-R2, 4-^4,1 ~ R2,lRiA 



R2,2Ri,i 


— i?2,4-R4,2 


R2,\Ri,2 


— i?2, 2-^4,1 


^2,2^4,4 


— i? 2 , 4-^4,2 


R2,s,Ria 


— i?2,2-R4 : 3 


-^2,2-^4, 4 


— R2,aRa,2 


-^2, 4-^4,3 


— R2,3,RiA 



R2,2Ra,A ~ R2,4R4,2 



(23) 



and 



N%(E) = 


L1.2L4.4 


— LiiL i2 


^2,2^4,4 


— L2. 4-^4,2 


A p L {E) = 


^3,2^4,4 


— -^3,4-^4,2 


L2, 2-^4,4 


— ^2,4^4,2 


N h L {E) = 


^2,2^3,4 


— L2, 4L32 


L2, 2-^4, 4 


— L2, 4-^4,2 


A h L {E) = 


L1.4L2.2 


— L12L2A 


-^2,2^4,4 


— ^2,4-^4, 2 



Using Eqs. (|23l24j) specified at E 
rewrite Eq. (f2"2")l as 



u 2 + 1 - {N P R N P L 



2ia.f 



+ A R A 



(24) 



0, we may then 



+ c.c.}w = , (25) 



with all the scattering amplitudes in Eq. ([25|) evaluated 
at E = 0. 

Putting together Eqs. (|21l25[) , one readily derives 
Eq.(9) and the related formula for 
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u +(x) 



u -(x) 



(21) 



with vp — ap/m and u±{x) being the roots of the 
second-degree equation 

+ GoAx)]u + G ltl = . (22) 
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